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A NOTE ON THE LADYZENSKAJA-BABUSKA-BREZZI CONDITION 



JOHNNY GUZMAN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS 



Abstract. The analysis of finite-eleinent-like Galerkin discretization techniques for the station- 
ary Stokes problem relies on the so-called LBB condition. In this work we discuss equivalent 
formulations of the LBB condition. 



1. Introduction 

The well known Ladyzenskaja-Babuska-Brezzi (LBB) condition is a particular instance of the so- 
called discrete inf-sup condition which is necessary and sufficient for the well-posedness of discrete 
saddle point problems arising from discretization via Galerkin methods. If X/j denotes the discrete 
velocity space and M/^ the discrete pressure space, then the LBB condition for the Stokes problem 
states that there is a constant c independent of the discretization parameter h such that 



^' (LBB) c||g,|U. < sup iilj^^, Vg, £ M,. 

The reader is referred to [5] for the basic theory on saddle point problems on Banach spaces and 
their numerical analysis. Simply put, this condition sets a structural restriction on the discrete 
^ ■ spaces so that the continuous level property that the divergence operator is closed and surjective, 

see [I] 3] , is preserved uniformly with respect to the discretization parameter. 

In the literature the following condition, which we shall denote the generalized LBB condition, 
is also assumed 



Cn! ^nTT^n^ IIV7 II ^ „ Ini'^-Vh)qh 

O 

(N 



(GLBB) c||V(?,,||l2 < sup ■"'. ' , ^Qh e M^, 



here and throughout we assume Alt C H^{fl). By properly defining a discrete gradient operator, the 
case of discontinuous pressure spaces can be analyzed with similar arguments to those that we shall 
present. Condition (jGLBB[) . for example, was used by Guermond ([Slin]) to show that approximate 
solutions to the three-dimensional Navier Stokes equations constructed using the Faedo-Galerkin 
j^ I method converge to a suitable, in the sense of Scheffer, weak solution. On the basis of (jGLBBp . the 

same author has also built ([10]) an H^'-approximation theory for the Stokes problem, < s < 1. 
Olshanskii, in [T2], under the assumption that the spaces satisfy (jGLBB[) carries out a multigrid 
analysis for the Stokes problem. Finally, Mardal et al., [TT], use a weighted inf-sup condition to 
analyze preconditioning techniques for singularly perturbed Stokes problems (see Section [S] below). 
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It is not difBcult to show that, on quasi-uniform meshes, (jGLBBp implies (jLBBI) . see [5]. We 
include the proof of this result below for completeness. The question that naturally arises is whether 
the converse holds. Recall that a well-known result of Fortin j^ shows that the inf-sup condition 
()LBB|) is equivalent to the existence of a so-called Fortin projection that is stable in Hj(ri). In this 
work, under the assumption that the mesh is shape regular and quasi-uniform, we will show that 
(IGLBBI) is equivalent to the existence of a Fortin projection that has L^-approximation properties. 
Moreover, when the domain is such that the solution to the Stokes problem possesses H^-regularity, 
we will prove that (|GLBBp is in fact equivalent to (ILBBp . again on quasi-uniform meshes. 

The work by Girault and Scott ([7]) must be mentioned when dealing with the construction 
of Fortin projection operators with L^-approximation properties. They have constructed such 
operators for many commonly used inf-sup stable spaces, one notable exception being the lowest 
order Taylor-Hood element in three dimensions. However, (jGLBBp has been shown to hold for 
the lowest order Taylor-Hood element directly [8]. Our results then can be applied to show that, 
(jGLBBp is satisfied by almost all inf-sup stable finite element spaces, regardless of the smoothness 
of the domain. 

This work is organized as follows. Section [2] introduces the notation and assumptions we shall 
work with. Condition (|GLBBp is discussed in Section |31 In Section 2] we actually show the equiv- 
alence of conditions (jLBBp and (jGLBBp . provided the domain is smooth enough. A weighted 
inf-sup condition related to uniform preconditioning of the time-dependent Stokes problem is pre- 
sented in Section [5l where we show that (|GLBBp implies it. Some concluding remarks are provided 
in Section [5] 

2. Preliminaries 

Throughout this work, we will denote by J7 C K'' with d = 2 or 3 an open bounded domain with 
Lipschitz boundary. If additional smoothness of the domain is needed, it will be specified explicitly. 
L^(f2), H^{n) and HqIH) denote, respectively, the usual Lebesgue and Sobolev spaces. We denote 
by L?g{n) the set of functions in L^(J7) with mean zero. Vector valued spaces will be denoted by 
bold characters. 

We introduce a conforming triangulation Th of 51 which we assume shape-regular and quasi- 
uniform in the sense of [2[. The size of the cells in the triangulation is characterized hy h > 0. We 
introduce finite dimensional spaces X/j C HQ(r2) and Mh C L?g{n) nH^{n) which are constructed, 
for instance using finite elements, on the triangulation 7/i. For these spaces, the inverse inequalities 

(2.1) \\vh\\H^<ch-^\\vhh2, ^Vh^y^h, 
and 

(2.2) hhWrn < ch-^hhh^, yqh 6 Mh, 

hold, see [2]. Here and in what follows we denote by c will a constant that is independent of h. 

We shall denote by Ch '■ Hj(i7) — > 'X.h the so-called Scott-Zhang interpolation operator ([13]) 
onto the velocity space and we recall that 

(2.3) \\v-Chv\\-L- + h\\Chv\\H^<ch\\v\\Hi, Vvenlin). 

and 

(2.4) \\v-Chv\\iii<ch\\v\\ti2, Vi.eHj(0)nH2(0) 

The Scott-Zhang interpolation operator onto the pressure space Xh '■ L?g{ft) — > Mh can be defined 
analogously and satisfies similar stability and approximation properties. We shall denote by nji : 
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L^(r2) — > X/i the L^-projcction onto X^ and by Ho : L'^{U) — > L'^{n) the L^-projection operator 
onto the space of piecewise constant functions, i.e., 



°°'%5m(X')*^' "'^"'< 



'{SI). 

I / I \ /rn / 

Ten 
For one result below we shall require full H^-regularity of the solution to the Stokes problem: 

Assumption 1. The domain il is such that for any f G L^(r2), the solution {ip,d) G Hj(il) X 
L? (j(i7) to the Stokes problem 

V-ip = 0, in n, 

^i = 0, on dV,, 



(2.5) 



satisfies the following estimate: 

(2.6) mu-^ + mm<4fh-- 

Assumption [T] is known to hold in two and three dimensions {d ~ 2, 3) whenever J7 is convex or 
of class C^'^, see [5J Theorem 6.3]. 

By suitably defining a discrete gradient operator acting on the pressure space, the proofs for 
discontinuous pressure spaces can be carried out with similar arguments. 

We introduce the definition of a Fortin projection. 

Definition 2.7. An operator J'h ■ Hj(i7) — > X/i is called a Fortin projection if Tf^ = Th and 

(2.8) / V-(« - Fhv)qH =0, Vv e Hj(fi), yqu e Ah. 

We shall be interested in Fortin projections T^ that satisfy the condition: 

(FHi) II-^/.«1Ihi<c1K>||hi, yvenlin), 

or 

(FL2) \\v~J'hv\\-L2<ch\\v\\iii, VveHlin). 

Let us remark that the approximation property (jFL2p implies H^-stability. 

Lemma 2.9. If an operator J^h : Hj(51) — > X;; satisfies (jFL2p then it is H^-stable, i.e., (jFHip is 
satisfied. 



Proof. The proof relies on the stability and approximation properties ()2.3p of the Scott-Zhang 
operator and on the inverse estimate (12.11) . for if u e Hq(J7), 

||-^/iw||hi < \\J^hV -ChvWiii +c||w||hi < ch^^WThV - Chv\\i^2 +c||w1|hi 

< ch^^\\v ~ Thv\\-L2 +ch^'^\\v ~Chv\\-L2 + c||u||hi- 

Conclude using the L^-approximation properties of the operators J-^ and Ch. O 

Remark 2.10. Girault and Scott, [7], explicitly constructed a Fortin projection that satisfies (JFH1|) 
and (jFL2p for many commonly used spaces. In fact, they showed that the approximation is local, 
i.e., 

\\J^hV - w||l2(t) + hrllThV - w||hi(t) < c/it||w||hi(a^(t)), Vu e Hj(f2) and VT G %, 
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where J\f{T) is a patch containing T. In particular, they have shown the existence of this projection 
for the Taylor-Hood elements in two dimensions. In three dimensions they proved this result for all 
the Taylor-Hood elements except the lowest order case. 

In this work we shall prove the implications 

(iLBBl) ^=^ 3Th s.t. (EH) and (IFHTI) 



(|GLBB|) <^^> 3Th s.t. dlH) and (|FL2|) •^= (iLBBj) and Assumption H 

thus showing that, in our setting, all these conditions are indeed equivalent. The top equivalence 
is well-known, see [HIHIIS]. The left implication is also known (see [5]), for completeness we show 
this in Theorem 13.31 The bottom implications, although simple to prove, seem to be new. 

3. The Generalized LBB Condition 

Let us begin by noticing that the generalized LBB condition (jGLBBp is actually a statement 
about coercivity of the L^-projection on gradients of functions in the pressure space. Namely, 
(jGLBBp is equivalent to 

(3.1) hh'^qhh^ > c||Vg,,||L2, Wqh 6 Mh. 

It is well known that (jGLBBp implies (jLBBI) . For completeness we present the proof. We begin 
with a perturbation result. 

Lemma 3.2. There exists a constant c independent of h such that, for all g/j G M^, the following 
holds: 

c\\qh\\L^< sup h^:^!^ + h\\Vqh\\^.. 



Proof. The proof relies on the properties (j2.3|) of the Scott-Zhang interpolation operator Ch, 

c\\qh\\L^- < sup II I, < sup ,, — TT] h sup -— — 

< sup — — — 7 1- sup —— , 

«,.ex,. ||Vu/,||l2 «GHi(n) IIVwIIl^ 

conclude using (|2.3p . D 

On the basis of Lemma [X^ we can readily show that (jGLBBp implies (jLBB[) . Again, this result 
is not new and we only include the proof for completeness. 

Theorem 3.3. (jGLBBp implies (ILBB]) . 

Proof. Since we assumed that Af/,. C i? ,,($7) fl H^{n), the proof is straightforward: 

sup -— 7 = sup — > -— = — 11 > C/I TTfeVgfe L2 

where, in the last step, we used the inverse inequality (j2.ll) . This, in conjunction with Lemma 13.21 
and the characterization p.ip . implies the result. D 
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Let us now show that the generahzed LBB condition (|GLBBp is equivalent to the existence of a 
Fortin operator satisfying (IFL2p . We begin with a modification of a classical result. 

Lemma 3.4. For all p e H^{n.) there is w G HQ(ri) such that 

Vv =p-IIqp, v\dT^O yT eTh, 

and 

XTeTh I 

Proof. Let p e H^{n) and T £%■ Clearly, 

P-Hop^ 0. 



It 
A classical result (pQ [TU [HI H]) implies that there is a wt S Hj(r) with V-vt = P — Hop in T and 

(3.5) I1Vut||l2(t) < c\\P - ^op\\l^{t)- 

Given that the mesh is assumed to be shape regular, by mapping to the reference element it is seen 

that the constant in the last inequality does not depend onT £Th. 

Let V € Hg(Sl) be defined as v\t — vt for all T in Th- By construction, 

V-w = p — LIoP, a.e. in fl. 
Moreover, 

TeTh TeTh TeTh Ten 

The first equality is by definition; then we applied the Poincare-Friedrichs inequality (since v\t — 
Vt G Hg(r)); next we used the properties of the function vt and the approximation properties of 
the projector Ho. D 

With this result at hand we can prove the following. 

Theorem 3.6. If there exists a Fortin operator J^h that satisfies (|FL2p . then (jGLBBp holds. 

Proof. Let qh G Mh . Using the properties of the operator XIq and the local analogue of the inverse 
inequality p.2p . we get 

l|Vgh||L2 = Y ll^(9'i~no(7h)||L2(T) < Y j^Wih-'^oqhWhf^T) <j^hh-'noqh\\l2. 

Ten Ten ^ 

From Lemma l3.4l we know there exists v G Hj(51) with V-w = qh — Hoq/i and 

1|v||l2 <c/i2||Vg,,||L2, 
hence 

W^QhWh <^\\qh-noqh\\l2^^ f {Vv){qh-noqh) = ^ f {V-v)qh, 

where the last inequality follows from integration by parts over each T and using the fact that 
v\dT — (see Lemma [3^ . 

Using the existence of the operator J^h, 

\\y<i,\\h < ^ fiy-^Hv)qh < ( sup h^:^) ^wt^vw^.. 
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It remains to show that 

\\Fhv\\i,2 <c/i2||Vg,,||L2. 
For this purpose, we use the approximation property (jFL2|) and Lemma 13.41 

\\:Fhv\\i,2 < \\Thv - v\\i^2 + \\v\\i^2 < c/i||vw||l2 + ch^wyqhh^ < ch^wyqhh^, 

where the last inequahty holds because of (|3.5p . D 



The converse of Theorem 13.61 is given in the following. 
Theorem 3.7. // (jGLBB|l holds, then there exists a Fortin projector J-^ that satisfies (|FL2p . 
Proof. Let v G Hj(ri). Define [zh^ph) E X/i x M^ as the solution of 



Zh-Wh - / Ph^-Wh = / v-Wh, ^Wh G X,i, 



In Jn 

Notice that (jGLBB|l provides precisely necessary and sufficient conditions for this problem to have 
a unique solution. 

Define J^hV := z/i we claim that this is indeed a Fortin projection that satisfies (JFL2I) . By 
construction, (12. 8p holds (see the second equation in (13.81) ). To show that this is indeed a projection, 
assume that v ~ v^ G X^ in (j3.8p . setting Wh ~ Zh ^ Vh we readily obtain that 

\\zh -Vh\\j,2 = 0. 

It remains to show the approximation properties of this operator. We begin by noticing that 
(|GLBBP implies 

(3.9) c||V,.||..< sup IflpL< sup ^"V?^^'"' ^"---^-'-"^-' 

w^ex^ \\whU2 wuey^H \\whK2 

where we used p.8p . To obtain the approximation property (|FL2p we use the Scott-Zhang inter- 
polation operator Cu, 

W^hV - u||l2 = / {ChV - v)-{ThV -v) + {ThV ~ Chv)-(FhV - v) 
Jn Jn 

< \\ChV - u||l2 \\ThV - v||l2 + / {J^hV - Chv)-{ThV - v). 

Jn 

We bound the first term using the approximation property (|2.3p of Ch. To bound the second term 
we use problem (|3.8p with Wh — J-^hV — ChV, then 

{ThV - Chv)-{FhV -v)= phV-{FhV - Chv) = / ph^-{v - Chv) = Vph-{v ~ Chv), 

In Jn Jn Jn 

we conclude applying the Cauchy-Schwarz inequality and using (|3.9p . D 
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4. Smooth Domains 



Here we show that, provided (jLBBp holds and, moreover, the domain Q, is such that Assumption!!] 
is satisfied, then (jFL2| holds and hence (jGLBB|) holds as well. This is shown in the following. 

Theorem 4.1. Assume the domain fl is such that the solution to (|2.5p possesses 'H? -elliptic regu- 
larity, i.e., Assumption\^holds. Then (ILBB[) implies that there is a Fortin operator J^h that satisfies 
(lFL2]l . 

Proof. Let v £ Hq(51). Define {zh,Ph) G X/j x Mh as the solution to the discrete Stokes problem 



Vzh.Vwh - I pjy-Wh = / Vu:Vw/i, Vw^ G X/j, 

QhVzh == / QhVv, yqh e Mh, 

In Jn 

where, in (I4.2p . the colon is used to denote the tensor product of matrices. Notice that (|LBBp 
implies that this problem always has a unique solution. 

Set ThV := Zh- Proceeding as in the proof of Theorem 13. 71 we see that this is indeed a projection. 
Moreover, p.8|) holds by construction. It remains to show that (|FL2|) is satisfied. To this end, 
analogously to the proof of Theorem 13. 7[ we notice that (|LBB[) implies 

lb/i|JL2 < c\\SI{ThV - w)||l2. 

We now argue by duality. Let ip and (j) solve (|2.5p with / = J-hV — v. Assumption ()2.6|) then implies 



Jn 

= f V{FhV - v) : V(^ - Chtp) - f {0- IhO) ViFhV - v) 
Jn Jn 

V{FhV - v) : V(C,,^) - / [Ihe] V-{FhV - v) 
Jn 



Notice that since 1^6 e M/j, j^^{IhO) V-{FhV — v) ^ 0. Since V-ip = 0, using (|4.2p . the estimate for 
Ph, dH and (EH), 



/ V{FhV--v) -.^iChij)^ [ phV-iChi' - i') < ch\\v - Fhv\\iii\\v - Fhv\\-L2. 
Jn Jn 

A direct application of of (|2.4p . (|2.3p and (|2.6p allows us to obtain the following estimates: 

f (9- Ih0) ViFhV ~v)+ f V{FhV ~ v):V{^ - Ch^) < ch\\FhV - v\\i^2 \\v\\ui 
Jn Jn 

We conclude using a stability estimate for (|4.2p 

\\J^hV - vWl-^ < ch\\FhV - v\\ii^ < ch\\v\\ii^, 
which, given (jLBBI) . is uniform in h. D 
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5. The Weighted LBB condition 



In relation to the construction of uniform preconditioners for discretizations of the time dependent 
Stokes problem, Mardal, Schoberl and Winther, [11], consider the following inf~sup condition, 

(5.1) c\\qh\\m+e-.L- < sup iil^^^^i^, Vg^£il4. 

where 

Ikll W.-i^ = ,^i^f^, (ll^ill?.. + e"'lk2|li.) , 

and 

By constructing a Fortin projection operator that is L^-bounded they have showed, on quasi- 
uniform meshes, that the inf-sup condition (|5.ip holds for the lowest order Taylor-Hood element 
in two dimension. In addition, they proved the same result, on shape regular meshes, for the mini- 
element. Here, we show that (|5.ip holds if we assume (|GLBBI) . A simple consequence of this result 
is that, on quasi-uniform meshes, (|5.ip holds for any order Taylor-Hood elements in two and three 
dimensions. 

Theorem 5.2. Let 51 be star shaped with respect to hall. If the spaces "K-h and Mh are such that 
(jGLBBp is satisfied, then the inf-sup condition (j5.ll) holds with a constant that does not depend on 
e or h. 

Proof. We consider two cases: e > h and e < h. 

Given that the domain J7 is star shaped with respect to a ball, we can conclude (JT]) that the 
following continuous inf-sup condition holds, 

(5.3) c\\q\\H^+,-iL^ < sup ..y^'"' , Vg e Lloi^), 

with a constant c independent of e. 

We first assume that e > ft.. Using (|5.3p for q^ e M^ we have, 

/ qh Vi^ J^^qhyjThv) WJ'hvh^neH^ 
c\\qh\\m+e-^L^ < sup -f^ = sup 



- ™P ir^\ ^^P -in\ ' 

i-hGX,, ||W/i||L2neHi t,eHi(0) ll^llL^ncHi 

where we used that, since (jGLBBp holds. Theorem 13.71 shows that there exists a Fortin operator Th 
that satisfies (|2.8p . By Lemma [^7^ and the approximation properties (jFL2p of the Fortin operator, 

\\J'hv\\-L2n,iii < c(|| J",,i;||l2 + e\\J^hv\\ii^) < c{\\v\\i,2 + \\v - Thv\\-L2 + e\\v\\-tii) 

< c{\\v\\i^2 + (e + h)\\v\\H^) < c{\\v\\i^2 + 2e\\v\\Hi) < c||w||L2neHi, 

where we used that h < e. 

On the other hand, if e < ft, we use gi = qh and 172 = in the definition of the weighted norm for 
the pressure space. Condition (jGLBBp then implies 

II II ^ ||V7 II ^ /n qh y-vh J^ qh V-f /» Ik/tlk^rK^Hi 
MhllH^+.-^L^ <c\\Vqh\\v^ <c sup -^ — <c sup --^ sup —r — . 



By the inverse inequality (|2.ip . 
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IVHU^neU^ <ca + £/»-^) 



Conclude using that e < h. D 

6. Concluding Remarks 

There seems to be one main drawback to our methods of proof. Namely, all our results rely 
heavily on the fact that we have a quasi-uniform mesh. However, at the present moment wc do not 
know whether this condition can be removed. Finally, it will be interesting to see if (|LBBp is in 
fact equivalent to (|GLBBI) on domains that do not satisfy the regularity assumption (|2.6p (e.g. non 
convex polyhedral domains). 

On the other hand, it seems to us that condition (jGLBBp must be regarded as the most important 
one. Our results show that, under the sole assumption that the mesh is quasi-uniform, this condition 
implies the classical condition (|LBB|) (Theorem 13. 3p . Moreover, as shown in Theorem 15.21 this 
condition implies the weighted inf-sup condition (|5.ip on quasi-uniform meshes. 
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